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Abstract 

This paper is devoted to the evaluation of the pionic scalar density at finite tem- 
perature and baryonic density. We express the latter effect in terms of the nuclear 
response evaluated in the random phase approximation. We discuss the density and 
temperature evolution of the pionic density which governs the quark condensate 
evolution. Numerical evaluations are performed. 



Introduction 



Pions play a crucial role in chiral symmetry restoration, due to their Goldstone boson character. 
The amount of restoration is measured by the modification of the order parameter, i.e. the quark 
condensate density, with respect to the vacuum value. In this context the evaluation of the expec- 
tation value of the squared pion field, linked to the scalar density of pions pj by = m 7r ($ 2 ), is 
of a great interest. For a single nucleon this quantity governs the total amount of chiral symmetry 
restoration of pionic origin, according to [1,2]: 

^ J d 3 x(N\1> 2 (x)\N) = 2m q J d 3 x(N\£*qq(x)\N) = £^ = (1) 



where N w is the scalar number of pions in the nucleon cloud, is the part of the nucleon S 
commutator of pionic origin and (N\A 7T qq(x)\N) represents the corresponding modification of the 
quark condensate with respect to the vacuum value. Similarly, in a uniform nuclear medium of 
density p, or in a heat bath, the evolution of the quark condensate originating from the pions is 
linked, to one-pion loop order, to the average value ($ 2 ) by : 

A'(gg(p,T)) (<fr 2 ) 
(OT(0,0)> 2/2 ' 1 ) 
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The same quantity ($ 2 ) governs also the quenching factor, 1 — ($ 2 )/3/J, of coupling constants 
such as the nucleonic axial coupling constant qa or the pion decay one f n , originating from pion 
loops, which is the counterpart of the mixing of axial and vector currents [3,4]. This type of 
quenching has to be seen as a manifestation of chiral symmetry restoration and should also apply 
to the case of p meson excitation by virtual photons as enters in relativistic heavy ion collisions. 
It is therefore interesting to evaluate the quantity ($ 2 ) in the conditions of such experiments. 
Now, the fireball which is the source of the dileptons contains, besides thermal pions, a significant 
residual baryonic background We have therefore to understand how the pion density evolves at 
finite values of both temperature and baryonic chemical potential. The first order approximation 
for the quantity ($ 2 ) adds the values for a pure heat bath and for a cold baryonic medium, writing 
with obvious notations : 

($ 2 )(p,T) = ($ 2 Mp = 0) + ($ 2 ) p (T = 0). (3) 



However this approximation is likely to be crude since the temperature has an effect on the pion 
density of the nuclear medium and on the other hand the presence of the baryonic background 
modifies the number of pions thermally excited. As an illustration of the second point, the pion 

density in the baryonic vacuum is fixed by the Bose-Einstein factor which is {e Wk ^ T — l) for 

pions of momentum k, with io k = \Jk 2 + rn^. In the nuclear medium the pion becomes a quasi- 
particle with a broad width. It can decay for instance into a particle-hole pair which has a smaller 
energy than the free pion. Its excitation is then favored by the thermal factor. There is therefore 
a mutual influence between temperature and density that we will investigate. 

The article is organized as follows. In the first section we evaluate the pion density of a nuclear 
medium at zero temperature. We relate this quantity to the nuclear response to a pion-like excita- 
tion. We evaluate this response in the RPA scheme, taking also into account the two particle-two 
hole {2p-2h) excitations. Within this framework, we study the deviation with respect to the inde- 
pendent nucleon approximation. In the second section we introduce the effect of the temperature 
through the modification of the nuclear responses. In the third section we incorporate the influence 
of the finite baryonic chemical potential in the heat bath case. 



1 Pion scalar density in the cold nuclear medium 



This quantity was discussed in relation to the quark condensate modification by Chanfray and 
Ericson [1]. They discussed its deviations with respect to free nucleons, introducing the nuclear 
response to pion-like excitation, treated in the static case — > oo. The extension to the non- 
static situation can be performed through the time-ordered graphs of Fig. 1 where the cross 
represents the point at which the pions are created or annihilated. The pion momentum is denoted 
q and oo is the excitation energy of the nuclear system in the intermediate state. The sum of the 
four graphs leads to the following expression which has also been derived with a different method 
in Ref. [5] : 
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where Q 2 = uo 2 — q 2 , v(Q 2 ) is the form factor of the pion vertex for which we use a monopole 
form : 
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Finally, Rl represents the spin-isospin longitudinal response function : 

R L (u, q) = £ | (n\ £ a t .q if |0) | 2 % - 25 n ). 
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Fig. 1. Time-ordered graphs for pion exchange. 

In the actual calculation, we add in the response Rl the excitation of the Delta resonance with 
the standard replacement : 



ff » </ T i — ► bl.q (T- )' 

fi , 7TiVJV 



(7) 



where 5 f (T f ) is the spin (isospin) transition operator connecting the spin (isospin) | and | states 
[6]. Note that for a matter of convenience we have incorporated in this operator the ratio of the 
ttNA and ttNN coupling constants. We have also assumed the same form factors, v(Q 2 ), at the 
7rNN vertex and at the 7riVA vertices. Finally we recall the link between the response function 
and the polarization propagator II (u;, q, q') which we will use in the following : 



V 

R(u,q) = hnU(u,q, q) 

71 



(8) 



We will first discuss the result for the free nucleon. 



1.1 Free nucleon. 



In the nucleon case, where the response Rl reduces to a simple expression, Eq. (4) provides for 
the (scalar) pion number N n in the nucleon cloud : 
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In the above equation we have used the following notations : e q = q 2 /2Mn and u& = Ma — 
M N + q 2 /2MA- The energy dependence of the Delta width Ta is taken from the analysis of the 
pion-nucleon scattering [6]. 

Our numerical inputs for the evaluation of N n are defined as follows. Without form-factor the 
integrals diverge linearly The resulting value is then quite sensitive to the cut-off function. We 
stress, however, that the specific case of the single nucleon is not the purpose of this paper. It 
is for us an element of comparison to introduce the nuclear effects. Since we want to restrict 
the calculation to a region where nuclear effects are reasonably under control, we have limited 
the integrals to : q — 1 GeV and ou — 1 GeV. The parameters are chosen so as to obtain a 
value = ^m^N^ = 30 MeV which is well in the accepted range [7-9]. This is achieved with 
A = 1 GeV and (g wNA /g„ NN ) 2 = 3.8. 



1.2 Infinite nuclear matter. 



We now turn to the case of infinite nuclear matter. In order to evaluate the response functions, we 
use the method of Delorme and Guichon [10] who calculated the zero order response in the local 
density approximation and then solved exactly the RPA equations. Their zero order response func- 
tion also includes the 2p-2h excitations. The corresponding Feynman diagrams of such processes 
are displayed in Fig. 2. 



Fig. 2. Feynman diagrams for the 2p-2h processes. The double line represents the Delta resonance. 

They are calculated in two steps. First we single out the contributions which reduce to a medium 
modification of the A self-energy, for which the parametrization of Ref. [11] is used. This parametriza- 
tion includes some 3p-3h excitation states as well. For the rest we use the results of Shimizu- 
Faessler [12] who evaluated the two nucleons p-wave pion absorption at threshold (cu = m n ), 
from which the A self-energy part, already taken into account, is separated out. As for each 
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the remaining contributions, an energy extrapolation suggested by the many-body diagrammatic 
interpretation is performed. 

We solve the RPA equations in the ring approximation : II = II + II V II. Here V is the particle- 
hole (p-h) interaction with the standard formulation : V = V n + V g > where the second piece is the 
short-range Landau-Migdal part. More explicitely with our definition of the response, V reads : 

v «W) / — <f — ,\ 

The corresponding Landau-Migdal parameters g' are different in the various channels : g 1 ' NN for 
the NN sector, g' AA for the A one, g' NA for the mixing of NN and AiV excitations. We adopt 
the following values: 

g' NN = 0.7, g' NA = g' AA = 0.5 . (11) 

The results are illustrated on Fig. 3 which shows the energy dependence of the zero order and 
RPA responses, for a fixed value of the momentum q = 300 MeV. The bare response presents a 
low-energy peak corresponding to the NN' 1 excitation and a high-energy one (AiV -1 excitations). 
The figure displays the RPA enhancement of the low energy peak, introduced by Alberico et al. 
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Fig. 3. Response function per nucleon at normal nuclear density as a function of the energy for a fixed 
momentum q = 300 MeV. The dot-dashed and continuous lines represent respectively the zero order and 
RPA responses. 

[13], which arises from the attractive nature of the p-h interaction. It also displays the collective 
behaviour of the A excitation with a splitting into two branches. 
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Once the energy and momentum integrations are performed, we find, as in Ref. [1], a moderate 
increase as compared to the free nucleon value, with the values per nucleon at normal nuclear 
density : 

t w N = Y? A jA = 38.5 MeV and N\j A = 0.55 (12) 



versus 30 MeV and 0.40 respectively for the free nucleon. The quantity is the (medium modi- 
fied) effective sigma commutator. 

More precisely, we decompose the response function II (accordingly S 71 ") into four types, depending 
on the kind of states which are excited at each external vertices : NN, NA, AN and A A (Fig. 4). 
The density evolution of the different components of the sigma commutator are represented on 




Fig. 4. Symbolic representation of the NN, AN, N A and AA response functions. The Delta resonance 
is represented by the double line. 

Fig. 5 both without and with RPA. Notice that, in absence of the RPA, there is already a contri- 
bution of the NA channel at finite density due to the 2p-2h excitations. The overall RPA increase 
of the sigma commutator mainly comes from that of the Ena + San parts. These last quantities 
embody the mixing of the AiV -1 configurations into the NN -1 ones. This is well known to be 
responsible for the enhancement of the low energy response (i.e. the NN~ X excitations) [14]. 



2 Inclusion of temperature 



2. 1 Influence on the virtual pion cloud 



We now introduce the temperature, via the Matsubara formalism : all the integrals over energy 
are replaced by an infinite sum over Matsubara frequencies. For the NN" 1 sector, the generalized 
Lindhart function (the imaginary part of which is proportional to the response function) at finite 
temperature is discussed in textbooks (see e.g. Ref. [15]). The result is the replacement of the 
Heaviside functions that characterize the occupation of fermion states at zero temperature by the 
Fermi-Dirac distributions. We have generalized this procedure to particle with width and applied 
it to A-N" 1 , N-A" 1 and A-A~ x rings. The generalized Lindhart function L ab -i for a process 
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Fig. 5. Density evolution of the different components of the effective sigma commutator: left (right) figure 
is without (with) RPA. 

involving a particle of type (a) and a hole of type (b) is found to be : 



L ab -i(ou,q) = — — - / k aka(cosv) r ■ . , — - — 



^ui + u) + ^{uj% +q -uj) 



fH){i-m +q )) 



(13) 



(14) 



where iV^/ is a constant arising from the summation of spin and isospin and T a ' b (uj) represents 
the width of the particle of type (a,b) for an energy uj. The occupation number of hadron species 
a is : 



/M = 



exp(^jA) + 1 



(15) 



where \i is the (common) chemical potential for baryons, the value of which fixes the baryonic 
density p at a given temperature. As implicitely stated before, we limit ourselves to nucleons and 
deltas i.e. p = piv + Pa- 

Concerning the A width, things are somewhat more complicated. In the medium, the pionic decay 
channel A — > n N is partly suppressed due to Pauli blocking. At the same time, other channels 
open, the pion being replaced by lp-lh, 2p-2h, etc... At normal density and zero temperature 
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the pionic channel remains dominant according to Ref. [11]. It represents approximatively 75% 
of the total width (~ 90 MeV to be compared with 120 MeV at the resonance energy) and the 
non-pionic decay channels only the remaining 25%. In view of the difficulties of a full calculation 
of the temperature effects, we have adopted the simplified following strategy : we have kept 
the parametrization of Ref. [11], derived at T = 0, for the non-pionic decay channel. We have 
introduced the temperature effects via the Matsubara formalism only for the main pionic part of 
the width. 

The Figure 6 displays the temperature evolutiom of the zero order response function. The dashed 
line, which represents the T = case, exhibits the N-N^ 1 and A-iV _1 structures. Increasing 
temperature tends to wash out more and more these peaks. At the same time an overall suppression 
effect occurs. Note that the A branch is less affected by the temperature because of the higher 
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Fig. 6. Zero order response at a baryonic density equal half nuclear matter density (p = 0.25m^.) as a 
function of energy for a fixed momentum q = 300 MeV for three temperatures. Dashed lines : T = 0. 
Thin full line : T = 50 MeV. Thick full line : T = 150 MeV. 

energies involved. In the RPA case (see Fig. 7) , one observes a similar behavior : an important 
general decrease and the loss of the lower energy structures. 

We now turn to the question of the thermal pions present in the heat bath. 
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Fig. 7. Same as before but for the RPA response. 
2.2 Inclusion of thermal pions. 

The effects previously discussed concerned only the modification, due to the temperature, of the 
virtual pion density present in the nuclear medium. At finite temperature, thermally excited pionic 
modes (quasi-pions) are also present and they give an additional contribution to the pion scalar 
density. Here for a better illustration we give the result for ($ 2 )/2/^ which according to Eq.(2) 
governs the amount of chiral symmetry restoration of pionic origin: 

($ 2 > p 3 gl NS f d 3 q J ( 1 1 \ 2 {n 2, o ( x 

ptl < $ 2 > T 

where n(u) = 1/ (exp(u/T) — 1) is the Bose occupation factor and D(u, q) the quasi-pion propa- 
gator. The second identity defines the quantities and < $ 2 >t, associated respectively with the 
first and second pieces of the r.h.s. of the above equation: represents an effective, temperature 
dependent, sigma commmutator per baryon, whereas < $ 2 >t is the scalar density of quasi-pions 
thermally excited. 
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3 Results 



In order to display the condensate evolution, all the forthcoming figures show its relative decrease, 
i.e. the quantity $ 2 /2/^ according to Eq. (2). We stress again that the points of interest are 
the influence of temperature on the nuclear pionic cloud contribution, the influence of the bary- 
onic density on the thermal pions one and finally how large is the deviation from the additive 
approximation of Eq. (3). 

We first present on Fig. 8 the contribution of the virtual pion cloud alone alone (term in 
Each box shows the temperature evolution at fixed baryonic density. It illustrates the suppression 
effect of the temperature which originates in the quenching of the nuclear response previously 
mentioned. The increase with baryonic density observed in this figure reflects the obvious fact 
that the pionic density follows the baryon one. 




Fig. 8. Relative decrease of the condensate coming from the virtual pion cloud part as a function of the 
density and the temperature. Each box corresponds to a fixed density as indicated. The density increases 
by steps of0.2po from left to right between 0.4 and 1.6po- In each box the points correspond to temperature 
increases by steps ofO.lOm^ between 0.05 and 1.05m,,-. 

In Figure 9 we present in the same fashion the condensate decrease due to thermal pions alone 
(term in ($ 2 )^/2/ 2 ). The iso-temperature curves show the influence of the baryonic density which 
pushes down the quasi-pion excitation energy, thus increasing their thermal excitation. 

Finally in the last figure (Fig. 10) we present the sum of both contributions. The competition 
between the variations of both terms with respect to the temperature is the source of the observed 
parabolic type shape. For comparison, we have plotted in open circles the approximation of Eq. (3) 
where the effect of the thermal pions at zero density is simply added to that of the pion cloud at 
zero temperature. This approximation does not display the hollow shape of the exact calculation: 
it overestimates the latter by at most 15% in the medium part of the temperature range we 
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Fig. 9. Same as the preceding figure but for the thermal pions. As a guidance in order to display the 
influence of the density on the thermal pions, points of equal temperature are joined by a dashed line. 



have considered, the deviation becoming quite small beyond T m 90-100 MeV. In this region the 
decrease of the pionic cloud contribution with temperature (Fig. 8) practically compensates the 
enhancement of the thermal excitations by density effects. The most important conclusion which 
can be drawn from Fig. 10 is that, due to the nuclear pions, the pion scalar density is much larger 
than in the absence of nuclear effects, already for densities of the order of 0.6 po- 




Fig. 10. Sum of both contributions of the virtual pion cloud and the thermal pions (black circles). For 
comparison the open circles represent the additive approximation of Eq. (3) 
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Conclusion. 



In conclusion we have studied the evolution of the quark condenaste of pionic origin under the 
simultaneous influence of the baryonic density and temperature. It is related to the scalar pionic 
density which comes on the one hand from the virtual nuclear pions and on the other hand from the 
thermally excited ones. We have expressed the first contribution in terms of the nuclear response 
to a pion-like excitation and evaluated it for the case of nuclear matter in the RPA scheme, 
first at zero and then at finite temperature. We have shown that the RPA produces a sizeable 
enhancement (~ 30%), while instead the temperature washes out the peaks and suppresses the 
nuclear response, hence decreasing the virtual pion density. 

As for the thermally excited pions we have shown that the presence of the baryonic background 
appreciably enhances their number. The cause has to be found in the lowering of the quasi-pion 
excitation energies, which favours their thermal excitation. When the densities of both types of 
pions are added, the mutual influences which go in opposite directions cancel their effects to a large 
extent. In the density and temperature domain that we have explored, the additive assumption of 
Eq. (3) which neglects the mutual influence is a good approximation. It deviates from the exact 
result by no more than 15%, the deviation being maximum around T m 50 MeV. At this T value 
the additive approximation slightly overestimates the pionic density. 

Our study has shown that, even at moderate baryonic density, the virtual nuclear pions are a 
major component of the overall scalar pion density. As an example, at nuclear matter density, 
they dominate in the temperature range we have considered, i.e. up to at least T as 150 MeV. 
Since the pion is the agent for the mixing of the vector and axial correlators, the consequence of 
our study is that the existence of a baryonic background, if any, should not be ignored in this 
mixing. 
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